Monopoles, confinement and charge localization in the t-J model with dilute holes 
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We present a quantum field theoretic description on the t-J model on a square lattice with dilute 
holes (i.e. near half-filling), based on the compact mutual Chern-Simons gauge theory. We show 
that, due to the presence of non-perturbative monopole plasma configuration from the antiferro- 
magnetic background, holons (carrying electric charge) are linearly confined and strongly localized 
even without extrinsic disorder taken into account. Accordingly, the translation symmetry is spon- 
taneously broken at ground state. Such an exotic localization is distinct from Anderson localization 
and essentially rooted in intrinsic Mott physics of the t-J model. Finally, a finite-temperature phase 
diagram is proposed. The metal-insulator transition observed in in-plane resistivity measurement 
is identified to a confinement-deconfinement transition from the perspective of gauge theory. The 
transition is characterized by the order parameter "Polyakov-line" . 

PACS numbers: 74.40.Kb,74.72.-h 



I. INTRODUCTION 

A major issue of the single-band t-J model on a square 
latticep] [2] is how the doped holes (i.e. Zhang-Rice 
singlets [3 modeling the copper-oxygen hybridization) in- 
teract with spin background eventually producing plenti- 
ful phase diagram of cuprate high-temperature supercon- 
ductors. Analytic study on this model is intricate due to 
the projective Hilbert space (i.e. two electrons are pro- 
hibited to simultaneously occupy the same site) in which 
electrons are "fractionalized" with spin-charge separation 
and emergent gauge degrees of freedom. [2j 0H6] These ex- 
otic phenomena open a new window for condensed matter 
physicists searching for unconventional quantum states of 
matter, and also shed lights on the amazing unified de- 
scription of fundamental physical laws at different length 
scales. [7] 

At half-filling (each site is occupied by one elec- 
tron), the t-J model is reduced to quantum Heisenberg 
model which possesses antifcrromagnetic long range or- 
der (AFLRO) and Mott insulating property. [1 [3] At 
this parent state, electric charge of electrons is totally 
frozen out rendering a pure spin model. Near half-filling 
where hole concentration is extremely low, charge local- 
ization and insulating ground state have been read out 
either directly or indirectly from various types of exper- 
imental justification. |10ffT3"] We know that disorder in- 
duces Anderson localization [TJ]. It is thus fundamen- 
tal to ask whether the localization near half-filling of 
cuprates is due to extrinsic disorder effect [10] or intrinsic 
quantum effect of the single-band t-J model[inj. In other 
words, is the translation symmetry at the ground state 
of the t-J model spontaneously broken? 

Toward this direction, new progress has been made 
recently. Z. Zhu et al. perform a large-scale density 
matrix renormalization group (DMRG) numerical sim- 
ulation by keeping enough number of states in each 
DMRG block with high accuracy. [16] The numerical re- 
sult demonstrates that the density of single hole is lo- 



calized in ladder systems, which was explained by the 
so-called phase string theoretical argument (More details 
and analysis can be found in Ref. I15j) . Experimentally, 
C. Ye et al. [T7] are able to significantly enlarge bias 
range in scanning tunneling microscopy (STM) to study 
the atomic scale electronic structure of the Ca2CuC>2Cl2 
parent compound with electron-donated defects. A well- 
defined in-gap state appears near the edge of upper Hub- 
bard band (UHB) and is strongly localized in which the 
typical localization length is order of lattice constant. [T7] 
These new findings combined with previous experimental 
hints call for a coherent theoretical description on charge 
dynamics of the t-J model near half-filling which is the 
main purpose of this work. 

In this paper, we shall provide a quantum field- 
theoretic approach to show that charge localization is 
driven by the non-perturbative monopole effect of a 
compact U(l) gauge dynamics that was initially in- 
troduced in Polyakov's seminal papers [IS]. We shall 
further propose the finite temperature phase diagram 
(Fig. [T]-a) near half-filling in which the metal-insulator 
transition is identified to monopole-driven confinement- 
deconfinement transition characterized by the order pa- 
rameter "Polyalov-line" jTSHST] . Essentially, electrons are 
fractionalized into bosonic holons (carrying charge) and 
bosonic spinons (carrying spin), both of which are mutu- 
ally entangled via a mutual Chern-Simons action. [221 [23"] 
Spinons are condensed in Arovas-Auerbach formalism [9 
and form superfluid at extremely low doping at zero tem- 
perature. At finite temperatures we assume spinons are 
still condensed (by adding a weak interlayer AF coupling 
J± or in-plane anisotropy a xy in spin space [8] \W \ I24H26] ) 
with a nonzero Neel magnetic transition temperature T/v 
where Bose condensation occurs, (see Fig. Ilia) In such a 
saddle point ansatz, we study charge dynamics at ground 
state and finite temperature behavior below T/y. 

First, the effective theory of charge dynamics is a 
(2+l)D compact U(l) gauge theory coupled to holons 
and ±27r phase-vortices arising from spinon superfluid. 
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introduced. The charge localization shall be derived in 
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FIG. 1. (Color online) (a) Intrinsic phase diagram near half- 
filling in absence of extrinsic disorder. T and 8 stand for 
temperature and hole concentration, respectively. The phase 
region of T < Tjv supports antiferromagnetic long range or- 
der (AFLRO). In this spin-ordered region, there is an addi- 
tional temperature scale Tmi about which the system under- 
goes metal-insulator transition. The high (low)-T region is 
metallic phase ( "Monopole-Plasma-Insulating Phase") as to 
be explained in main texts. The dashed line segments im- 
ply that the present analysis is in fact self-consistent only 
near half-filling, (b) A schematic description of holon local- 
ization in AF background. Black arrows form AF spin back- 
ground and blue balls are two holons. The red directional 
circle stands for a "-27r" phase-vortex whose core is occupied 
by a holon. The three objects inside the dashed circle form a 
localized bound state. 



Since the gauge dynamics is in confined phase where 
monopole operators are relevant, holons and ±2tt phase- 
vortices are enforced to form gauge-neutral bound states. 
Because phase-vortices are static particles with infinite 
effective mass, the bound states eventually localize even 
without extrinsic disorder taken into account. Hence, the 
translation symmetry is indeed broken spontaneously. 

Second, a confinement-deconfinement transition of 
charge degree of freedom occurs at a finite temperature 
Tmi- At T < Tmi, holons are confined via linear potential 
and localized such that charge transport is of insulating 
nature. This low-temperature phase is called "Monopole- 
Plasma-Insulating Phase" as to be explained below. At 
T > T"mi j the linear confinement disappears, but holons 
still perceive a logarithmic interaction from static phase- 
vortices of spinon superfluid background. Therefore, this 
high-temperature metallic phase supports a metallic be- 
havior. Tmi is identified to the metal-insulator transition 
temperature scale observed in in-plane resistivity mea- 
surement of heavily underdoped cuprates. [101 |2"T1 12"5] The 
phase diagram is shown in Fig. [l]-a. Disorder is absent in 
the present theoretical approach but physically disorder 
effect in a realistic material is expected to further amplify 
localization and thereby enhance Tmi- 

The paper is organized as follows. In Sec. [TTJ the 
compact mutual Chern-Simons gauge theory of the t-J 
model is reviewed where frequently used notations are 



Sec. Ill and the discussion of finite temperature phase di- 
agram is arranged in Sec. |IV| All results are summarized 
in Sec. N\ 



II. COMPACT MUTUAL CHERN-SIMONS 
THEORY OF DOPED MOTT INSULATORS 

In this section, we first review our understanding of 
doped Mott insulators, especially the phase string effect 
of the t-J model and its compact mutual-Chern-Simons 
gauge theory. More details can be found in Refs. |2"51 
[29ti32l We stress that the present t-J model only contains 
nearest-hopping term and indeed the electron- and hole- 
doped cases are symmetric although the realistic cuprates 
have asymmetric phase diagram. By considering more 
hopping terms and super exchange terms, it is possible 
that the asymmetric phenomena can be realized which is 
beyond our present theoretical analysis. 

Let's start with the so-called "sign structure" of doped 
spin model. The significance of figuring out the sign 
structure of a given theory can be illustrated by the Na- 
gaoka state [33 , one of the few exact results about the 
t-J model. The Nagaoka problem, the U = oo Hubbard 
model with one hole, is equivalent to the large J limit 
of Kondo lattice model with one conduction electron[34]. 
It can be shown that the validity of Nagaoka theorem 
relies heavily on the triviality of the sign structure of 
one hole hopping term of the t-J model^T. The sign 
structure of a model is called trivial when the Hamilto- 
nian matrix has only non-positive off-diagonal elements 
in some well-chosen basis. If a Hamiltonian has trivial 
sign structure and satisfies some other conditions, we can 
use Perron-Frobenius theorem to show that the ground 
state of the system is non-degenerate and has a positive 
wave function in this well-chosen basis. It is the triv- 
ial sign structures of the AF Heisenberg model and the 
one-hole infinity-[/ Hubbard model that give rise to the 
Marshall theorem and the Nagaoka theorem. 

At finite doping and finite J, however, the sign struc- 
ture of the t-J model is highly nontrivial[35] . We can also 
choose the Marshall bases {|<^>)} to make the off-diagonal 
matrix elements of Heisenberg term non-positive. Nev- 
ertheless, nontrivial signs appear in the hopping term, 
resulting in a model where the Perron-Frobenius theo- 
rem is no longer applicable. This can be seen easily in 
slave-fermion formalism, where the hopping and Heisen- 
berg terms can be written as 30J 

H t = -t <Tflf j b] a b ia + h.c, (1) 



(2) 



Here, bi a is the annihilation operator of bosonic spinon 
at site i with spin a (a — ± denote spin up/down), and 
fi the annihilation operator of fermionic holon at site 
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i. The configurations with more than one electrons at 
a certain site are projected out, leaving a Hilbert space 
with constraint 

E 6 Uv+4/i = i- (3) 

a 

The electron operator is written in this formalism as 

Cia = (-<y) l f}b lai (4) 

where the Marshall sign (— a) 1 depends on the sublattice 
index (i even or odd) and the spin a of the electron. The 
sign structure of the Heisenberg term Eq. ^ is trivial, 
while the spin dependent sign structure of hopping term 
Eq. ([I]) indicates that a minus sign appears whenever a 
down spinon exchanges with a holon. This is called the 
phase string effect[30\. 

To keep track of this this nonrepairable phase string 
effect, we define a nonlocal unitary transformation 

e^^exp^X^COnNuj, (5) 

where 6i(l) — Imln(zi — z;), and zi = xi + iyi is the 
complex coordinate of site I. Under this nonlocal unitary 
transformation O —> e ze Oe~ , the t-J model becomes 

H t = -t 2 e iA h:h,,- TV //„/;,„ • h.C, (6) 

hj = - j - £ j-MA-^'^-^- (7) 

{i,j),o,o' 

We have also used the Jordan- Wigner transformation in 
2D to describe holons by bosonic operators hi's instead of 
fermionic operators fiS. The sign structure of the origi- 
nal t-J model after the nonlocal unitary transformation 
is now represented by two compact U(l) gauge fields A s 
and A h defined by 

A h = £ E ft(0 - 9 M « - O ( mod 2 *)> (8) 

A% = \Y,M l )- e Mv>l (mod27r). (9) 

Those equations indicate that the holons (spinons) are 
the 7r (±7r) vortices of the gauge field A h (A s ). 

We can now use mean field theory to deal with this 
model, since the phase string effect is explicitly tracked 
by the nonlocal unitary transformation Eq. ([5]). This 
procedure bears a resemblance to the simple direct prod- 
uct variational wave function treatment of the Haldane 
phase of ID AF spin-1 chain after a nonlocal unitary 
transformation, which makes it possible to use local fer- 
romagnetic order parameter to reveal the nonlocal hid- 
den string order of the phsse[3S]. The validity of the 
above method can be also illustrated by the ID t-J model 



which possesses non-Fermi-liquid behavior. In ID case, 
9i(l) — it ■ 9(1 — i), where 9(x) is the Heaviside step func- 
tion. According to the definitions Eq. Q and (|9|, the 
gauge fields Af^ and A\^ vanish. Thus the two terms in 
the t-J model Eq. ^ and Eq. ^ both have a triv- 
ial sign structure, with all the signs absorbed into the 
definition of the fractionalization of physical electron op- 
erator. The simple mean field treatment is then enough 
to get the correct Luttinger-liquid behavior of correlation 
functions [30], 

At mean field level, the t-J model is reduced to the 
effective phase string model, with Hamiltonian H e ff = 
H h + H s , and 

H h = -t h ^ h\ h 3 + h - c -» ( 10 ) 

{i,j) 

H s = -J s £ e^&l&J_ ff +h.c, (11) 

where, th and J s will be defined in below. This model 
possesses a compact U(l)<g>U(l) gauge symmetry: 

fci->M M S (12) 

b ia ^b ia e ia9 >, (13) 
A s,h^ A s,h + d s,h_ e s,h^ (u) 

It is instructive to formulate the above mean field the- 
ory in path integral formalism after regularizing appro- 
priately on a cubic spacetime lattice. This is the so-called 
mutual Chern-Simons gauge theory 22, 23 . Each holon 
(spinon), as a n (±7r)-ffux tube, constitutes "electromag- 
netic" flux of A 1 ^ (A^) which is minimally coupled to 
spinon (holon). This structure leads to compact mutual 
Chern-Simons topological term in which the compact- 
ness of A s ^ h is exactly protected. In such a formalism, 
the partition function encompasses compact U(1)<8>U(1) 
gauge redundancy: 

Z= J2 I ' D[A s ,A h ]D[h\h,b\b]e~ s (15) 

in which the action S — ^ x C (x denotes spacetime co- 
ordinates) with C = Ch + C s + £mcs- £-n (£ s ) is the La- 
grangian density for holon (spinon) , which includes usual 
gauge covariant operators responsible for minimal cou- 
pling between holon (spinon) and A s (A h ). For further 
derivation, we need to explicitly write down the concrete 
form of C s : 

C s =bl (d - zaA>> + X s ) b irT + ^ (b\ a b la ) 2 

-J s (e™<b\ +&<J b\_ a + h.c) , (16) 

where, Einstein summation is employed for all indices, 
a, /?,••• denote space directions x or y. The indices 
/i, v, A • • • stand for spacetime directions. d v is forward 
difference operator on 3-dimensional spacetime lattice. 



4 



U2 is onsite repulsion energy which softens hard-core con- 
dition of Schwinger bosons. X s is chemical potential. 
The effective AF superexchange energy J s is defined as: 
Jg = t;A s where bosonic resonating- valence bond (RVB) 

order parameter A s = Y^a( eiaA ™b\ +&a b\_ (J ) ^ below 
the pseudogap temperature T [37]. This meanficld RVB 
condensate gaps out the usual U(l) gauge fluctuation 
(denoted as XJ(l) A ) which arises from single-occupancy 
constraint, in contrast to the strong gauge fluctuation in 
U(1) A slave-boson theory[5j. In the latter, it is necessary 
to carefully investigate confinement-deconfinement of the 
U(1) A gauge fluctuation by integrating out high energy 
modes of matter fields 2, 38 . The compactified mutual 
Chern-Simons term is expressed by: 

£mcs = ~ ^ K - d„ (A h x - 2*J$) , (17) 

where, is antisymmetric tensor of rank-3. ■yV' t f' h are 
two integer-valued link variables to take account of the 
periodicity of the gauge field A s ' h . 



III. MONOPOLES, CONFINEMENT AND 
LOCALIZATION 

A. Emergent compact gauge theory 

In this section, we focus on the AF phase in the for- 
malism of mutual Chern-Simons theory. A compact U(l) 
gauge field theory emerges in the limit of low hole con- 
centration. 

At half-filling, spinon condensation leads to AFLRO. 
Near half-filling, we still assume spinon condensation, i.e. 
(ba) 7^ 0. Let's formally write spinon field b a = {^/n^ + 
h)e where Uq is (^-dependent condensation fraction. By 
integrating out the massive amplitude fluctuation h of 
spinon field [39], one can obtain the following effective 
Lagrangian: 

+ V" A (AJ - 2^)d v {A\ - 2«JI>?) + C h , (18) 

where, g = J s /^, gi = .92 = ^n J s . jV£ plays the 
key role of the static ±2ir phase- vortices arising from the 
spinon superfluid. 9 is a scalar function of spacetime 
coordinate. By further employing the unitary gauge |4"0], 
9 can be absorbed into A^, while, £mcs keeps invariant 
due to the antisymmetry property of . Since d^9 £ R, 
one obtains a massive real vector field A h u with A'^ £ M. 

Therefore, the first two terms in Eq. ( p~8| ) are replaced 
by " f (A%) 2 + 3 f(A^) 2n by keeping A^£ R in mind. 
Gaussian integration over A^ leads to: 

£=^(F° v -2nn tw f+£ c , (19) 



where, e = y/4noJ s is coupling constant of "emergent 
(2+l)D compact gauge dynamics". The "speed of light" 
is implicit here without loss of generality. The gauge field 
strength tensor F^ v is defined as F^ v = d^A s v — d v A s ^ (d M 
is backward difference operator on 3-dimcnsional space- 
time lattice). The new plaquette variable n^ u is defined 
as: rv = {d^Jf* -d v Jf*). C c = -i2e^ x A^Jf* + L h . 
The definition of here is locally well-defined. As will 
be discussed below, n M „ induces monopole configurations 
of three-dimensional Euclidean spacetime when jV^ has 
singularity. It is particularly interesting that in the quan- 
tum paramagnetic phase of non-linear sigma model, com- 
pact gauge degree of freedom is also found and plays an 
important role in classifying quantum spin liquids. [4"lTt44] 
In these systems, the compact gauge degree of freedom 
is essentially due to the Wess-Zumino-Witten topological 
term of quantum SU(2) spins. [U] 

It is clear that two kinds of particles simultaneously 
carry the gauge-charge of A s , namely, holon h in C h 
and static ±2tt phase-vortex arising from spinon super- 
fluid. Let's define (2+l)D current of phase-vortices as: 
<* vX d v jY£ = where, £ Z, J a = (static vor- 
tices). As such, C c can be simplified to 

C c = -iAt{2J Q + h)h) + C h [A$=0] , (20) 

where, [a s =o\ stands for Lh without Aq term. Hence, 
we find that the gauge-charge of ^fo is ±2 while each 
holon carries +1 gauge-charge, such that, only possible 
negative gauge-charge comes from 

B. Monopole plasma configuration 

Monopole effect generically gets suppressed by finite 
density of matter field. Here, the holon matter field which 
couples to A s is extremely dilute, such that the monopole 
effect is expected to be relevant. [43] Let's briefly follow 
Polyakov's approach[18 by using the present mathemat- 
ical symbols and explicitly keep track of monopole ef- 
fect in context of doped antifcrromagnets. In analog 
to the point-like solution "Dirac monopole" in three- 
dimensional space, we can define the "magnetic field" 
vector as below (in dual lattice): SS^ = \e tJ ' uX n v \ ■ 2ir. 
The divergence of 38 ^ is in general quantized at 2ir, i.e., 
dn3$n = 2irq where q is integer scalar field defined on 
spacetime sites, q 7^ if JV^ has singularity. In general, 
can be globally defined and factorized into longitu- 
dinal and transversal components separately j!8j (in the 
continuum limit, it becomes the Hodge decomposition of 
a general differential form into exact, co-exact and har- 
monic forms on a Riemannian manifold): 

n^u = [d^rriu + Xv) - d v (m^ + x M )] - e^ uX d\<f>, (21) 

where, m v is integer vector field, x„ is a real vector field 
with absolute value smaller than 1. <f> is a real scalar field. 
It is a linear equation on lattice and one can check that 
the degrees of freedom of those fields on both sides of the 
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equation are indeed the same. Substituting factorization 
formula of n MJ/ into ( 19 1 we find that the can be com- 
bined with m M + Xfi rendering A s „ G R. This rearrange- 
ment brings convenience for the Gaussian integration for 
A s . The effective Lagrangian is thus transformed into: 



2tt 2 

-—-fyd^d^ + C c , 



4g2^ i» 

where, A^ € R. In deriving above expression, the cross- 
ing term u F^ v t ,LvX dx4>" is neglected, because 



F, 



e^ x d x cb =2e^ x d^A s J^ = 2e^%(A s J^ 



where the first term is trivial boundary term and the 
second term vanishes due to antisymmetry property of 
e^ x . To proceed further, let's substitute factorization 
formula of n M „ into d^SS^ resulting in: — = q where 
A = g? m c£ m is lattice Laplacian. The formal solution of <f> 
can be written as: <fi — — A -1 */ . 

Finally, we obtain the final effective action: 



4e 2 



2tt 2 



X 

(22) 



where, the second term in the above action describes 
a 3D plasma of monopoles with Coulomb interaction 
(— A -1 )^' ~ \x — x'\ . The configuration {q} repre- 
sents distribution of point-like "magnetic charge" (i.e. 
monopole). This Coulomb gas (monopole plasma) rep- 
resentation of compact U(l) gauge theory serves as the 
starting point of our following discussions. 



C. Confinement and localization at zero 
temperature 

The monopole plasma has far-reaching consequences: 
the long-range interaction in monopole plasma spoils out 
the correlation of original gapless A^-photon at weak- 
coupling limit and generates a gap in the low-lying charge 
excitation spectrum. This gap generation can be viewed 
as an alternative physical picture of Mott physics which 
strongly freezes charge degree of freedom near half-filling. 

The matter field h can be neglected at half-filling. In 
the absence of C c , Eq. (19) is identified to (2+l)D pure 
compact U(l) gauge theory [15) 135]. with Eq. (22) as its 



monopole plasma representation. To show the confining 
nature of the (2+l)D pure compact U(l) gauge theory 
at zero temperature, we can introduce a scalar field \ 
to reexpress the instanton part of the partition function. 
Taking into account only q x — ±1 configurations, we 
obtain the following sine-Gordon action: 



S e s = 



]T((Vx) 2 -M 2 cosx) 



(23) 



where M 2 = (2ir/e) 2 exp(— const. /e 2 ). x pl avs the role of 
scalar potential of the Coulomb charges, and its gradient 
is the electric (magnetic) field. The appearance of small 
mass M of \ in weak coupling e limit leads to a short 
range correlation function of electromagnetic field of the 
original U(l) theory. It is the finite density monopoles 
with long range interaction that spoil the correlations. 

In order to probe confinement of gauge-charge, one can 
define Wilson loopjTT] as 



W[C] 



where, C is a temporal rectangular r x t loop with r 
(t) spatial (temporal) distance. is unit current and 
forms the directional loop C. The underlying potential 
between test-particle and test-antiparticle is defined as 
V(r) = — Hindoo l n ^ [ c l _ Repeating the same transfor- 
mation from Eq. (19) to Eq. (22), we find the similar 
formula for W[C]: 

W[C] r.Z-^exp ( 



x exp 



2?ri 



Qx ^)x,x' Qx 



(24) 



It describe monopole plasma with a fixed external 
monopole configuration Q x . The screening of the ex- 
ternal monopole configuration requires a free energy pro- 
portional to the area of the rectangle C. In other words, 
the Wilson loop exhibits area law at large t limit: W[C] ~ 
e~ Ktr with positive coefficient k for any given coupling 
constant e. Since strong coupling limit is always a con- 
finement state, we can reasonably draw the conclusion 
that (2+l)D pure compact U(l) theory is always con- 
fined at zero temperature. As a results, all virtual par- 
ticles that carry gauge-charge in the vacuum of (2+l)D 
pure compact U(l) gauge theory must be confined into 
gauge-charge neutral bound state. 

Near half-filling, the doped dilute holes may be directly 
viewed as test-particles in the gedanken-experiment 
which physically interprets Wilson loops jUj: resulting 
in linear confinement between holes and appropriate 
amounts of phase- vortices of spinon superfluid. The 
infinite effective mass of the latter leads to strong lo- 
calization of holes (carrying charge degree of freedom). 
As shown in Fig[l]-b, two holes form a localized state 
whose wave function may be expressed as |two holes) = 
• , • ) , where, -0: and • denote — 27r phase-vortex 
and hole, respectively. It carries zero gauge-charge (- 
2+1+1=0). As spin degree of freedom is energetically 
expelled away from phase- vortex cores, a holon must oc- 
cupy a core under the single-occupancy constraint of the 
t-J model, justifying Fig. [T]-b. Overall, we obtain the 
strong localization without extrinsic disorder/impurities 
taken into account. Existence of monopole effect in a 
spinful charge-neutral superfluid drives the charge local- 
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ization, and, the translation symmetry is broken sponta- 
neously. 

Since the present quantum field-theoretic approach is 
based on the phase string decomposition of electrons [301 
HH] as shown in Sec. |nj we emphasize that the essen- 
tial origin of localization mechanism can be traced back 
to the singular phase string effect discovered by Weng 
et al. [HI [301 EE] ■ The phase string effect is mathemat- 
ically captured by the exact "sign structure" of the t-J 
model by formulating partition function with the world- 
line path-integral. [35] Pictorially, the worldlines of one 
holon and one spin-! spinon wrap each other can con- 
tribute a minus sign under Gutzwillcr projection (the to- 
tal particle number of spinons and holons at each site 
must always be one). Such a particular sign structure 
encodes underlying non-Fermi liquid behaviors and im- 
plies the notion of "Sign Matter" introduced by Zaanen 
and Overbosch.l49l 



We conclude that finite temperature effectively reduce 
one dimension of our original theory. The Coulomb gas 
has only one plasma phase in 3D; while in 2D, there 
is a Berezinskii-Kosterlitz-Thouless transition from dis- 
order phase to critical phase. The sine-Gordon repre- 
sentation Eq. (23 1 of the Coulomb gas is also valid. 



And the 2D sine-Gordon model with inverse tempera- 
ture /?' = l/(/?e 2 ) is a well studied model in conformal 
field theory. At low temperature (large (3, small f}') of 
our original U(l) theory, the monopole effect term cos% 
is relevant, resulting in a short range correlation function 
of x, and the compact U(l) theory is in the confinement 
phase. At high temperature (small (3, large /?'), however, 
the monopole term is irrelevant, and the compact U(l) 
theory is in the deconfinement phase. The confinement- 
deconfinement transition do have measurable effect in our 
original t-J model. 



IV. PHASE DIAGRAM OF FINITE 
TEMPERATURE 

A. The effect of finite temperature 

At finite temperature, the compact U(l) theory has 
a deconfinement phase where the behavior of holons 
changes dramatically. We identify this confinement- 
deconfinement transition with the metal-insulator tran- 
sition observed in experiments. 

The same procedure dealing with zero temperature 
U(l) theory is valid at finite temperature. The two dif- 
ferences are: (i) Three dimensional infinite lattice is re- 
placed by a lattice with imaginary time size (3 = 1/kgT. 
And only periodic configurations along this direction con- 
tribute to the partition function, (ii) When calculat- 
ing the Green's function G(x — x') = — (A -1 )^/ in 
Eq. ( [22| , the integral along the temporal direction is re- 
placed by the Matsubara frequency summation. Using 
the fact that the Fourier transformation of the lattice 
Laplacian operator A x , x > = Y^^x,x'+n 
is 4^ sin 2 (fc / _ 1 /2), the Green's function in real space is 
given by 



G(r,x) 



i r d 2 k 



Eexp(zk • x + iu n T) 
- , , T 4V.sin 2 ^+4sin ; 



,2 cjj, 
2 



(25) 

One can show that at zero temperature limit (3 — > oo, 
the Green's function reproduces the three dimensional 
Coulomb potential at large distance 



G(r,x) 



f 



Vx 2 + 



= , (3^oo. 



(26) 



On the other hand, the high temperature limit f3 — > 
give us an two dimensional Coulomb potential 



G(x) 



1 r d 2 k exp(ik-x) 

(2tt) 2 4£ i sin 2 (fc l /2) 



, p^O. (27) 



B. Metal-insulator transition and order parameter 

Experimentally, magneto-resistance measurement [50 
of heavily underdoped YBa 2 Cu306+a; compounds indi- 
cates that the development of AF order has little ef- 
fect on the in-plane resistivity. It is also found that the 
magnitude of in-plane resistivity is so large that conven- 
tional band theory breaks down.[2Hl 02] In the present 
theory, gauge-charge is linearly confined at all e's at 
zero temperature. According to Svetitsky-Yaffe univer- 
sality arguments [55], it is suggested that a confinement- 
deconfinement transition exists at a finite temperature 
Tmi- Therefore, it is naturally addressed that electric 
charge transport may be explicitly altered when the sys- 
tem undergoes the transition. At T < Tmi holons are still 
linearly confined and strongly localized as same as ground 
state, implying insulating nature of charge transport. 
The only source for providing charge mobility is thermal 
fluctuation. At T > Tmi, the linear confinement disap- 
pears such that the charge mobility is enhanced. Albeit 
disappearance of linear confinement, the logarithmic in- 
teraction now plays the leading role and renders a metal- 
lic behavior. It is widely believed that the confinement- 
deconfinement transition temperature Tmi as a function 
of e starts from origin (0, 0) in T-e plane and roughly 
monotonically increases with the increase of e. |52l[55] 
Doping holes in general depletes the spinon condensa- 
tion fraction such that no(5) decreases with the increase 
of doping (5, so does the coupling constant e which is 
defined as y/4noJ s . Correspondingly, Tmi should mono- 
tonically decrease. Consequently, we identify Tmi as the 
metal-insulator transition temperature scale observed in 
in-plane resistivity measurement. [28 The order parame- 
ter of this transition is so-called "Polyakov-line" [TM2T) : 



O 



exp 



A° dr 



(28) 
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with = ksT and fee the Boltzmann constant. O — 
(^)O if T < (>)Tmi, which characterizes the preserva- 
tion (breaking) of central group. The Polyakov-line was 
ever utilized in "Short-Range Order phase" of Hubbard 
model by Wiegmann[56 but in quite a different context. 
Finally, we obtain the phase diagram shown in Fig. [T]-a. 
Starting from the present understanding on localization, 
quantitative study of charge transport is quite interesting 
and will be extensively addressed in our future work. 



C. Dual-type "Nernst" effect: A new quantum 
phenomenon 

At the end of the discussion, we in particular empha- 
size the novel vortex-core structure near half-hlling. In 
Fig. [l]R, a quantized unit electric charge (i.e. holon, the 
blue ball) is surrounded by spinful supercurrent. For the 
lower holon, the supercurrent is counter-clockwise; for 
the upper holon, the existence of the extra —2tt phase- 
vortex (red directional circle) leads to net clockwise su- 
percurrent. In analog to Nernst effect [57] in which vor- 
tex is formed by electric supercurrent, we predict that 
there is a dual-type effect if one can polarize spinful 
vorticity along z-direction (cuprate sample is put in xy- 
plane). Then, by applying temperature gradient along 
indirection, one can measure net spin accumulation at 
the two edges of y-direction. 



V. SUMMARY 

In conclusion, the present work provides a semi- 
quantitative quantum field-theoretic analysis on the long- 
standing problem: can electric charge be intrinsically lo- 
calized in the t-J model?. The monopole plasma config- 



uration which comes from the antiferromagnetic back- 
ground is proved to play fundamental role of driving 
charged holes into localized states. Although it is a tech- 
nical challenge to determine the concrete range of param- 
eter tj J in which the result makes sense from the present 
semi-quantitative analysis, our work has proved that the 
pure t-J model itself has the intrinsic ability to stabilize 
such translation symmetry breaking phase without the 
help of external disorder /impurity, which fundamentally 
differentiates the present localized ground state from An- 
derson localization. This result is a reasonable answer 
to the recent STM experimental finding and consistent 
to the phase string argument for DMRG numerical sim- 
ulation of ladder systems. Apart from this result, we 
also figure out the finite temperature phase diagram and 
especially theoretically explain the mechanism of metal- 
insulator transition found in electric in-plane resistivity 
measurement, which sheds light on a new way to reorga- 
nize transport experimental findings in curates in a single 
framework. A much more quantitative study along the 
present perspective is important and will be leaved to 
future work. 
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